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1 Introduction 



Little needs to be said nowadays about the importance of one-loop calculations. During the 
last years methods which (partially) automatize these calculations have become available |], 
|2], |H . The automatization is though nowhere as fully developed as for tree- level calculations 
0, [5]. The reason is that loop calculations are in general more involved, and moreover a 
complete automatization limits the possibilities to use a program beyond its designed scope. 
In virtually all implementations of one-loop calculations on the computer, dimensional 
regularization || [7]] has been employed for calculating renormalized amplitudes, and not 
without reason as it is a very elegant formalism that preserves most symmetries. However, 
dimensional regularization presents problems in chiral theories, related to the extension 
of 75 to arbitrary dimensions |8], |9|, |K], [11]]. In the Standard Model, a naive scheme — 
anticommuting 75 without further modifications — is usually used, but the situation is not 
completely clear ||12|| . More problematic are supersymmetric theories, where a variant, 



regularization by dimensional reduction [13], is commonly used in spite of possible in- 



consistencies at higher loops |14 |. On the other hand, a new method in 4 dimensions, 
constrained differential renormalization (CDR), has recently been developed [15|, I6| . It 



is a version of differential renormalization |L7j that preserves gauge invariance at least to 
one loop WE, [TJ], IS). There are some hints that it also preserves supersymmetry | 19fl . At 



any rate, it is clearly convenient to have at hand alternative methods implemented for 
automatized one-loop calculations. 

We have implemented the three mentioned methods in two program packages, called Form- 
Calc and LoopTools, for the evaluation of one-loop amplitudes. It turns out that CDR 
and dimensional reduction are equivalent at the one-loop level (at least in the cases where 
the latter is well-defined). This is a new result and will be discussed in the next sec- 
tion. Actually, once the CDR coordinate-space expressions and procedure are translated 
into momentum space, both methods can be set up in an identical way. Therefore, our 
programs have in fact only two different options: to perform calculations in dimensional 
regularization or in dimensional reduction / CDR. We shall refer to these two possibilities 
as calculating in D or 4 dimensions, respectively. 

Our programs do not fully automatize one-loop calculations to the extent that they directly 
produce cross-sections from the process specification. Instead, FormCalc produces output 
which can easily be evaluated further, e.g. numerically in a Fortran program. One big 
advantage of FormCalc is that it leaves the user with a Mathematica expression which is 
considerably easier to modify than Fortran code. LoopTools, on the other hand, supplies 
all the functions needed for the numerical evaluation of the FormCalc output. We supply 
also a demonstration program to show how this numerical evaluation can be done. 
The paper is organized as follows: In Section [| we outline the regularization methods we 
have used in the programs, and discuss the equivalence between dimensional reduction and 
CDR. In Section the implementation of the program packages is explained, where Section 
37TI describes the analytical and Section [3.2| the numerical part. In Section |] we illustrate 
for the example of elastic Z Z scattering the individual steps of a one-loop calculation using 
our packages. Section § lists the computer requirements needed to compile and run the 
packages, and their availability. Finally, an appendix collects the functions implemented 
for numerical evaluation. 
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2 Calculations in 4 and D dimensions 



Dimensional regularization has proved to be the most convenient method for comput- 
ing quantum corrections in gauge theories. The definition of the regularized expressions 
has two parts: analytic continuation of momenta (and other four-vectors) in the number 
of dimensions, D, and an extension to D dimensions of the Lorentz covariants (7^, g^, 
etc.). The second part is achieved by treating the covariants as formal objects obeying 
certain algebraic identities 0. The dimensionally regularized Feynman graphs are mero- 
morphic functions of the complex parameter D, and the poles at D = 4 can consistently 
be subtracted to all orders in perturbation theory. Such a minimal subtraction scheme 
(MS) defines a renormalization fulfilling the usual requirements of a quantum field theory: 
causality, unitarity, field equations, Ward identities, etc. R |TD[. Problems only appear for 



identities that depend on the 4-dimensional nature of the objects involved. This is the 
case for the Fierz identities or for relations using the Levi-Civita tensor. In particular, 
the extension of 75 to D dimensions is problematic. The consistent prescriptions of Refs. 
B H [10[ lead to spurious anomalies that have to be corrected with finite counterterms. 



We follow instead Ref. |12] and work with an anticommuting 75. We also maintain the 
usual relations for traces, such as Tr(7 5 7^7 !/ 7 P 7 cr ) = e^ upa Trl. Even though they can be 
incompatible with an anticommuting 75 in D dimensions, the resulting ambiguities are 
expected to cancel in non-anomalous theories fi"2| , |20[| . For one- loop calculations one can 
proceed in the following manner (see Ref. [2T| for details): 



1. Calculate traces and simplify the Dirac algebra in D dimensions. 

2. Write everything in terms of scalar and tensor integrals in D dimensions. 

3. Decompose tensor integrals into Lorentz-covariant tensors constructed from the ex- 
ternal momenta and the metric tensor. One is left with scalar integrals as coefficients 
of these tensors. 

4. Calculate the scalar integrals and tensor coefficients in D dimensions, expand the 
whole diagram in e = D — 4, subtract the - poles and take the limit e — > 0. Computer 
algebraically it is convenient to do this in two steps: 

(a) Add local terms for products of D times a divergent integral: DI — > AI+c, with 
/ the integral and c the coefficient of its - pole . Then, make D = 4 outside the 
integrals. 

(b) Calculate the scalar integrals and tensor coefficients in D dimensions, expand 
them up to order e°, and subtract their poles. 

A modified version of dimensional regularization, designed to preserve supersymmetry and 



gauge invariance, was proposed by Siegel |R|. It is based on dimensional reduction 



from 4 to D dimensions: while the integration momenta are D- dimensional, as in usual 
dimensional regularization, all other tensors and spinors are kept 4-dimensional. Therefore 
one works with two kinds of objects: D-dimensional and 4-dimensional ones. For the 
validity of the field equations and gauge invariance, one must impose the identity 

9^9 up = 9;, (1) 
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where is the 4-dimensional metric = 4) and g^ the .D-dimensional one (g M M = D). 
The Dirac algebra, including 75, is performed in 4 dimensions. Unfortunately, regulariza- 
tion by dimensional reduction is known to be inconsistent Hl~4| |. Nevertheless the incon- 



sistencies arise at higher orders and the method has successfully been applied to many 
calculations in supersymmetric theories.^ At one loop, the procedure described for dimen- 
sional regularization can be followed. The only differences are that the Dirac algebra is 
4-dimensional and that a 4-dimensional metric tensor and Eq. (p]) have to be introduced. 
Notice that g^ u only appears in the decomposition of the tensor integrals. These (and the 
tensor coefficients) are identical in dimensional regularization and dimensional reduction. 
Finally, differential renormalization is a method of regularization and renormalization 
in coordinate space that cures UV divergences by substituting badly-behaved expressions 
by derivatives of well-behaved ones. The method has proved to be quite simple and conve- 
nient in a number of applications. A symmetric procedure of differential renormalization 
preserving gauge invariance at the one-loop level has recently been proposed. This so-called 
constrained differential renormalization proceeds in two steps: 

1. The expressions read from Feynman diagrams are written in terms of a complete set 
of (singular) basic functions. The basic functions are products of propagators with 
derivatives acting on one of them. Functions are treated differently depending on 
whether their indices are self-contracted or not. All algebra is performed strictly in 
4 dimensions. 

2. These basic functions are substituted by their renormalized expressions. The renor- 
malization of the singular basic functions has previously been fixed once and for all, 
such that a set of simple rules is respected (these rules are just natural extensions of 
mathematical identities among tempered distributions) [TB |. In particular, the rules 



imply that renormalization does not commute with index contraction. 

Although differential renormalization naturally works in coordinate space, it is possible to 
perform the reduction of Step 1 in momentum space and use the Fourier transforms of 
the renormalized basic functions in Step 2, which correspond to the tensor integrals in the 
dimensional methods. Expressions and manipulations in coordinate space can easily be 



translated into momentum space (see Appendix B of Ref . JT(| ) . 

As mentioned in the introduction, it turns out that CDR and dimensional reduction are 
equivalent methods for one-loop calculations. The equivalence can easily be understood 
once CDR has been translated into momentum space. First, the minimally subtracted D- 
dimensional tensor integrals are identical in the limit D — > 4 to the Fourier transforms of the 
corresponding renormalized basic functions of CDR (which we shall also call tensor integrals 
in the following), up to a redefinition of the renormalization scale.0 The reason is that one- 
loop integrals in D dimensions satisfy the relations imposed by the CDR rules. Hence, a 
possible discrepancy can only arise in the initial conditions (namely, the renormalized value 
of the scalar one- and two-point functions), but this can be taken care of by the freedom 



1 Dimensional reduction can also be used in non-supersymmetric theories, but more care is required in 
order not to break unitarity p3[ . 

2 We have explicitly checked that the CDR tensor integrals in Appendix B of Ref. Jl6| coincide with 
the tensor integrals in dimensional reduction in the MS scheme if the CDR renormalization scale M is 
redefined as logAf 2 = log/1 2 + 2 where [i is the renormalization scale in the dimensional method. 
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to choose the renormalization scale. Second, all the algebra outside tensor integrals is 
considered 4-dimensional in both schemes. Third, the distinction of tensor integrals with 
contracted and uncontracted indices in CDR is imitated in dimensional reduction by Eq. 
(0) and the addition of local terms in Step |4a]. To illustrate this, consider the renormalized 



tensor integral (see Ref. [|lj for its definition), which is the same in both methods (up 
to the freedom to choose the renormalization scheme). In CDR the tensor integral is 
considered independently; its renormalized value differs from g^ v times the renormalized 
Cfiv by a local term fl5j . In dimensional reduction we have 

g^C^ = g» v I g^Coo + PH*P3" C H ) = Voo + (PiPj) C v 

V ij=i / ij=i (2) 



= DC 00 + ^ (PiPj) C ij = 4C oo - g + (P&j) C V ■ 

i,j=l i,j=l 

Cqo and Cy coincide in both methods. The extra local term, —1/2, is precisely what is 



needed to obtain C^u from g^ v Cn V in CDR [15, 16). That the same occurs in the general 



case follows from the fact that in dimensional reduction one can contract the 4-dimensional 
metric with D-dimensional integration momenta before performing the integrals, and the 
resulting contracted tensor integrals also satisfy the CDR relations. Notice that conven- 
tional dimensional regularization introduces extra D's coming from g^^s outside the tensor 
integrals. Also, the Feynman rules in dimensional regularization contain D in some the- 
ories. Hence it can render different results, and not only in intermediate steps. A simple 
example where the results in D and 4 dimensions differ is the electron self-energy in QED: 

dim - eg - ( ~ [4m eJ B (A; 2 ,m e 2 ,0) + 2^ 1 (A; 2 ,m e 2 ,0) + ^-2m e ] , (3) 




167T 1 



CDR, 2 
dim. red. 6 



[4m eJ B (fc 2 ,m 2 ,0) + 2^ J B 1 (fc 2 ,m 2 ,0)] . (4) 



167T 2 

The functions B and B\ are defined in App. [A|. 

Since dimensional reduction and CDR are equivalent at the one-loop level, one only needs 
a single implementation for the two methods. We have decided to use the CDR idea and 
add local terms when there are self-contracted indices after all the simplifications have been 
performed. This allows to work completely in 4 dimensions. Of course, the dimensionally 
regularized scalar functions and tensor coefficients are also used in this case. 



3 Implementation of the program packages 

The evaluation of one-loop diagrams using our packages proceeds in two steps: 
In FormCalc the symbolic expressions for the diagrams as obtained from FeynArts [p4 
simplified algebraically such that the output can be used almost directly in a numerical 
program, e.g. in Fortran or C++. 

LoopTools is the second of the packages and supplies the actual numerical implementations 
of the one-loop functions needed for programs made from the FormCalc output. It is based 



arc 
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on the reliable package FF [25] and provides in addition to the scalar integrals of FF also 
the tensor coefficients in the conventions of Ref. pifl . LoopTools offers three interfaces: 
Fortran, C++, and Mathematica, so most programming tastes should be served. 



3.1 FormCalc 

FormCalc is a Matiiematica-based program to calculate one- loop Feynman diagrams in 
either CDR or dimensional regularization. The program reads input from FeynArts and 
returns output in a way well suited for further numerical (or analytical) evaluation. Form- 
Calc was designed to handle large numbers of diagrams, for example running the complete 
W + W~ — > W + W~ amplitudes (about 1000 diagrams |26f ) takes approximately 10 minutes 
on a Pentium PC. 

The structure of FormCalc is rather simple: it prepares the symbolic expressions of the 
diagrams in an input file for FORM, runs FORM, and retrieves the results. This interaction 
is transparent to the user. FormCalc combines the speed of FORM with the powerful 
instruction set of Mathematica and the latter greatly facilitates further processing of the 
results. The following diagram shows schematically how FormCalc interacts with FORM: 



FormCalc 






FORM 

(very fast, but 
not so user friendly) 


OneLoop [amp] 


- ReadForm — *- 


Mathematica 

(user friendly, slow) 







user interface internal FormCalc calls 

The main function in FormCalc is OneLoop (the name is not strictly correct since it works 
also with tree graphs). It is used like this: 

« FormCalc .m; 
$Dimension = 4; 
amps = << myamps.m; 
alldiags = OneLoop [amps] 

The file myamps.m is assumed here to contain amplitudes generated by FeynArts. The 
dimension — D for dimensional regularization or 4 for dimensional reduction / CDR — is set 
with $Dimension. 

Alternatively, if one wants to evaluate only a subset of diagrams, 

somediags = OneLoop [ Pick [amps, {3, 5, {21, 29}}] ] 
will calculate diagrams 3, 5, and 21 through 29, or for a single diagram it can be just 

onediag = OneLoop [ amps [ [9] ] ] 

Note that OneLoop needs no declarations of the kinematics of the underlying process; it 
uses the information FeynArts hands down. When amplitudes are loaded, FormCalc auto- 
matically puts the external particles on shell and for some common processes^] introduces 

3 Currently processes with 4 or 5 external legs, i.e. 2 — > 2, 2 — » 3, 1^3, 1 — » 4. 
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kinematical invariants (e.g. Mandelstam variables for a 2 — > 2 process). The former can 
be disabled by setting $0nShell = False before loading the amplitudes. 
For example, the QED electron self-energy in Eq. ([3D was calculated with the following 
short program: 

<< FormCalc .m; 
$0nShell = False; 
$Dimension = D; 

OneLoop[« electronSE. amp] //. Abbreviations [] 
with the result 

EL 2 ME B0[Pair[k[l] , k[l]] , ME2, 0] 



4 Pi 2 

EL 2 Bl[Pair[k[l] , k[l]] , ME2 , 0] ga[k[l]] _ EL 2 (-2 ME + ga[k[l]]) 

8 Pi 2 16 Pi 2 

Similarly, Eq. (f|) can be reproduced by putting $Dimension = 4. This results in 

_EL 2 ME B0[Pair[k[l] , k[l]] , ME2 , 0] _ 
4 Pi 2 

EL 2 Bl[Pair[k[l] , k[l]] , ME2 , 0] ga[k[l]] 

8 Pi 2 



Even more comprehensive than OneLoop, the function ProcessFile can process entire files. 
It splits the results up into a bosonic and a fermionic part; this is commonly needed for 
the numerical evaluation e.g. if one wants to sum over fermion generations. ProcessFile 
is invoked e.g. as 

ProcessFile ["vertex . amp" , "results/vertex"] 

It reads the FeynArts amplitudes from vertex, amp and produces re suits /vert ex. m which 
contains the bosonic part and results/vertexF .m which contains the fermionic part. 
OneLoop and ProcessFile return expressions where spinor chains, dot products of and 
Levi-Civita tensors contracted with polarization vectors have been collected and abbrevi- 
ated. A term in such an expression may look like 

C0i[ccl, MW2, S, MW2, MZ2 , MW2, MW2] * 

(P12*S*(-8*a2*MW2 + 4*a2*MW2*S2 - 28*a2*CW~2*MW2*S2 + 
16*a2*CW~2*S*S2 + 4*a2*C2*MW2*SW~2) + 
047*S*(-32*a2*CW~2*MW2*S2 + 8*a2*CW~2*S2*T + 8*a2*CW~2*S2*U) - 
P13*S*(-64*a2*CW~2*MW2*S2 + 16*a2*CW~2*S2*T + 16*a2*CW~2*S2*U) + 
03*(-2*a2*MW2*T + a2*MW2*S2*T - 7*a2*CW~2*MW2*S2*T + 
4*a2*CW~2*S*S2*T + a2*C2*MW2*SW~2*T + 2*a2*MW2*U - 
a2*MW2*S2*U + 7*a2*CW~2*MW2*S2*U - 4*a2*CW~2*S*S2*U - 
a2*C2*MW2*SW"2*U)) 
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The first line stands for the tensor coefficient function Ci(M^, s, M^, M|, M^, M^) (see 
Appendix) which is multiplied with a linear combination of abbreviations like 047 or P12 
with certain coefficients. 

These coefficients contain the Mandelstam variables S, T, and U and some short-hands for 
parameters of the Standard Model, e.g. a2 = a 2 or C4 = cos -4 ^- Replacing such com- 
monly appearing factors by symbols makes expressions faster to handle for Mathematica. 
The coefficients are collected in a way optimized for numerical evaluation: for example, 
Mathematica will automatically add together the terms in -2*a2*MW2*T + a2*MW2*S2*T - 
7*a2*CW~2*MW2*S2*T once the numerical values for the various parameters are substituted 
(i.e. the result is of the form (a number) *T), so it is not worth simplifying this expression 
further. If desired, FormCalc can nevertheless perform even these last simplifications to 
obtain the most compact analytical form.f] 

The abbreviations like 047 or P12 are introduced automatically and can significantly reduce 
the size of an amplitude. Consider P13 in the excerpt of code above. It is composed of 

P13 = Oil - 012 - 013 + 014 + 020 

The Onn are in turn made up of 

011 = e23*kl3*k43 

012 = el3*k23*k43 

013 = el3*k24*k43 

014 = el2*k31*k43 
020 = e23*kl4*k43 

and the enn and knn are defined as 

el2 = Pair[e[l] , e[2]] 
el3 = Pair[e[l] , e[3]] 
kl4 = Pair[e[l] , k[4]] 

Pair [a, 6] denotes the scalar product of two four- vectors, and e[z] and k[i] the polar- 
ization vectors and momenta, respectively. 

To get an idea of how advantageous the introduction of abbreviations can be, it is useful 
to compare the Leaf Count of the expressions in Mathematica. The leaf count gives a 
measure of the size of an expression, more precisely it counts the number of subexpressions 
or "leaves" on the expression tree. P13 has a leaf count of 1 since it is just a plain symbol. 
In comparison, its fully expanded contents have a leaf count of 113. 

The definitions of the abbreviations can be retrieved by Abbreviations [] which returns 
a list of rules such that 

result //. Abbreviations [] 

gives the full, unabbreviated expression. It is of course necessary to save the abbreviations 
before terminating a session, e.g. with 

Abbreviations [] >> abbr 

4 FormCalc wraps two functions, ol and o2, around coefficients at different levels in the amplitude. The 
default setting ol = o2 = Identity can be changed e.g. to ol = o2 = Simplify. 
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3.2 LoopTools 



LoopTools supplies the functions needed for the numerical evaluation of the code produced 
by FormCalc. LoopTools follows the conventions of [21 1. The complete reference to the 
LoopTools functions can be found in Appendix [A|. 

LoopTools consists of three parts: a Fortran library, a C++ library, and a MathLink 
executable that can be used with Mathematica directly. 



3.2.1 Using LoopTools with Fortran 

To use the LoopTools functions in a Fortran program, the two files tools . F and tools. h 

must be included, tools .F contains actual Fortran code and should be included only once 

per program, tools .h contains the declarations and common blocks and must be included 

in every function or subroutine in which the LoopTools functions are used. 

Before using any LoopTools function, the subroutine bcaini must be called. At the end 

of the calculation bcaexi may be called to obtain a summary of errors. 

A very simple program would for instance be 



#include "tools. F" 



program simple_program 
#include "tools. h" 



call bcaini 

print *, BOCIOOODO, 50D0, 80D0) 

call bcaexi 

end 



Note that, as for all preprocessor commands, the # must stand at the beginning of the line. 
Several default values can be superseded by defining preprocessor variables before including 
tools. F, e.g. the renormalization scale can be changed with #define MUDIM ... or the 
IR regulator mass can be changed with #def ine LAMBDA . . . ; some more technical options 
are described in the manual [p7fl. 



3.2.2 Using LoopTools with C++ 

To use the LoopTools functions in a C++ program, the file ctools.h must be included. 
Similar to the Fortran case, before making the first call to any LoopTools function, 
bcaini () must be called and at the end bcaexi () may be called to get a summary of 
errors. 
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In C++, the same simple program looks like 

#include <f stream. h> 
#include "ctools.h" 

main() 
{ 

bcaini () ; 

cout « B0C1000., 50., 80.) « endl; 
bcaexi () ; 



The renormalization scale and the IR regulator mass can be changed with the functions 
set_mudim and set_lambda, respectively. 

3.2.3 Using LoopTools with Mathematics 

The Mathematica interface is probably the simplest to use: 

In[l] := Install [ M bca M ] 

0ut[l]= LinkObject[bca, 1, 1] 

In [2] := BO [1000, 50, 80] 

Out [2]= -4.40593 + 2.70414 I 

One-loop functions containing non-numeric arguments (e.g. B0 [1000, MW2, MW2]) remain 
unevaluated. 

Again, the renormalization scale and the IR regulator mass may be changed with the 
functions Mudim and Lambda, respectively. 



4 An application to the Standard Model 



We demonstrate in this section the basic steps of a real one-loop calculation using our 
packages. For an example, we choose the process ZZ — > ZZ in the electroweak Standard 



} 



Model || . 



A one-loop calculation generally includes 



Step 



Program 



typical CPU time 



1. Generate Diagrams 

2. Simplify analytically 

3. Produce Fortran code 

4. Compile with driver program 



FeynArts 
FormCalc 
(NumPrep) 
(num . F) 



3 min. 



3 min. 



7 min. 



10 min. 
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The quoted execution times are for the full one-loop ZZ — > ZZ calculation (about 500 
diagrams) as run on a standard Pentium PC under Linux. 

The numerical evaluation of the FormCalc output is not fully automated. The reason is 
that at one-loop level already many features appear which must be treated differently from 
process to process. The two programs NumPrep and num.F are supplied with the demo 
code to show how the numerical evaluation can be done. Although both were designed for 
2 — * 2 gauge-boson scattering processes, it should not be too difficult to adapt them at 
least to other 2 — > 2 processes. (Basically, one has to supply the kinematics in num.F.) 
In the following we present only a brief but characteristic excerpt of code for each step. 
The complete demo code is contained in the online distribution of FormCalc. 

Step 1: Here the generation of the self-energy diagrams is demonstrated: 

« FeynArts.m 

tops = CreateTopologies [ 1, 2 -> 2, (* 1 loop, 2 -> 2 *) 

ExcludeTopologies -> (* create only self-energies *) 

{Tadpoles, WFCorrections , Triangles, AllBoxes} ] 

inss = InsertFields [ tops, 

{V[2], V[2]} -> {V[2], V[2]}, (* V[2] = Z *) 

Model -> "SM", InsertionLevel -> {Particles}, 
Restrictions -> {NoGeneration2, NoGeneration3} ] 

amps = CreateFeynAmp [inss] ; 
ToFAlConventions [amps] » zzzz . self . amp 

This code excerpt shows a typical application of FeynArts: create the topologies, 
insert fields into them (V[2] is a Z boson in FeynArts lingo), and create the ampli- 
tudes, i.e. apply the Feynman rules. The function ToFAlConventions converts the 
full FeynArts 2.2 format back into the simpler FeynArts 1 format which is needed by 
FormCalc. 

Step 2: The so-created self-energy diagrams are then simplified with FormCalc. Calcu- 
lating in D and 4 dimensions yields equivalent results in the case of ZZ — > ZZ, hence 
we have omitted the explicit definition of $Dimension here (the default is D). 

« FormCalc. m 

ProcessFile ["zzzz . self . amp" , "self"] ; 
Abbreviations [] >> abbr 

Steps 3 and 4: The FormCalc results need to be converted into a Fortran program. (The 
numerical evaluation could, in principle, be done in Mathematica directly, but this 
becomes very slow for large amplitudes.) The simplest way to do this in Mathematica 
is something like FortranForm [result] >> file . f . With our demo code we supply 
a much more sophisticated program called NumPrep which goes well beyond simple 
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translation into Fortran code. For example, it groups the one-loop integrals into 
angle-dependent and -independent integrals, so that the latter need to be calculated 
only once e.g. when integrating over the angle at a fixed energy. 

NumPrep is too complex to describe here in detail. Instead, we give an example of 
how the final Fortran code looks like: 

#include <defs.h> 

double complex function self() 

implicit logical (a-s,u-z) 

implicit double complex (t) 
#include <vars.h> 

self = reso**2*(-5.7477663296703327914881096651780*ab61*03 - 
3.696762828539980268636174052200372*ab62*03 + 
ab59* (14. 787051314159921074544696208801489 + 

0.00022228809903668380734374191468819973*MH2)*03 + 
ab58* (22. 991065318681331165952438660711951 + 

0.0004445761980733676146874838293763995*MH2)*03 + 
0.0006668642971100514220312257440645992*ab60*MH2*03 + 
. 0020005928913301542660936772321937976*abl5*MH2**2* 
03 - abl4*03* (-215176. 58454926418740557299134203950 - 

0.00022228809903668380734374191468819973*MH2**2 + 

7.393525657079960537272348104400745*S) - ... 

Admittedly, this code looks ugly to the human eye, but note that it is highly optimized 
for numerical evaluation: apart from the 03 which is one of the abbreviations already 
introduced by FormCalc, the one-loop integrals have been replaced by variables (e.g. 
abnn for A and B functions). Also, the possibly resonant Higgs propagator 1 / (s — Mg) 
has been replaced by the variable reso so that it can be treated more easily in Fortran. 
These are just examples of what one can do with the FormCalc output. 

Of course, NumPrep produces also the code to calculate the abbreviations and a 
Makefile to conveniently compile the code. The Fortran code produced by NumPrep 
needs in addition a driver program which supplies it with the necessary parameters, 
kinematics, etc. This driver program is called num.F and is included in the demo. 



5 Computer requirements and availability 

FormCalc and LoopTools should compile and run without change on any Unix-based plat- 
form. In particular, they have been tested on DEC Alpha, HP 9000, and Linux. FormCalc 
needs Mathematica 2.2 or above including the MathLink compiler (mcc) and FORM 2 or 
above.Q LoopTools requires a Fortran- 77 compiler, the GNU make utility, and the GNU C 
and C++ compilers (gcc, g++). 

The programs FormCalc and LoopTools can be obtained via WWW from 



http : //www-itp . physik . uni-karlsruhe . de/f ormcalc| and 



http : //www-itp . physik . uni-karlsruhe . de/looptool4 respectively. 



FormCalc runs also with FORM 1, but cannot fully simplify spinor chains with external fermions then. 
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The packages contain a comprehensive manual giving installation instructions and a de- 
tailed description of every function. 

FeynArts is available from 

|f tp : / /ftp . physik . uni-wuerzburg . de/pub/hep/ index . html . 
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A Reference of the LoopTools functions 
A.l One-point functions 



Function call 


Description 


AO(ms) 


one-point function 



The real argument ms (double precision in Fortran) is the mass squared: 



ms = 




A. 2 Two-point functions 



Function call 


Description 


B0(ps, mis, m2s) 
Bl(ps, mis, m2s) 
BOOCps, mis, m2s) 
BlKps, mis, m2s) 


scalar two-point function 
coefficient of p^ 
coefficient of g^ u 
coefficient of p^Pu 



All arguments are real (double precision in Fortran) and are related to the momenta 
and masses as follows: 

ps =p 2 
mis = ml 
m2s = m\ 
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A. 3 Derivatives of two-point functions 



Function call 


Description 


DB0(ps, mis, m2s) 
DBKps, mis, m2s) 
DB00(ps, mis, m2s) 
DBlKps, mis, m2s) 


derivative of BO 
derivative of Bl 
derivative of BOO 
derivative of Bll 



All derivatives are with respect to the momentum squared. The arguments are as in the 
case of the two-point functions. 



A. 4 Three-point functions 



Function call 


Description 


CO (pis, p2s, plp2s, mis, m2s, m3s) 
C0i(id, pis, p2s, plp2s, mis, m2s, m3s) 


scalar three-point function 
three-point tensor coefficients 



Except for the id all arguments are real (double precision in Fortran) and are related 
to the momenta and masses as follows: 

pis = pj 
p2s = p\ 
plp2s = (pi +p 2 ) 2 
mis = m\ 
m2s = m\ 
m3s = m\ 

COi is a generic function for all three-point tensor coefficients. A specific coefficient is 
selected with the id argument: 

COiCccO, . . .) =C (...) 
COiCccOO, ...) =Cbo(...) (8) 
C0i(ccll2, . . .) = C 112 (. . .) etc. 

Since the indices are symmetric, the identifiers are assumed to be ordered, i.e. there is only 
ccl22 but not cc212. 
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A. 5 Four-point functions 



Function call 


Description 


D0(pls, p2s, p3s, p4s, plp2s, p2p3s, 

mis, m2s, m3s, m4s) 
D0i(id, pis, p2s, p3s, p4s, plp2s, p2p3s, 

mis, m2s, m3s, m4s) 


scalar four-point function 
four-point tensor coefficients 



Except for the id all arguments are real (double precision in Fortran) and are related 
to the momenta and masses as follows: 

pis 
p2s 
p3s 
p4s 
plp2s 
p2p3s 
mis 
m2s 
m3s 
m3s 

DOi is a generic function for all four-point tensor coefficients. A specific coefficient is 
selected with the id argument: 

D0i(dd0, . . .) = £>„(. . . ) 
DOKddOO, . . .) = Ax)(...) (10) 
D0i(ddl223, . . .) = L>i223(. • • ) etc. 

Again, since the indices are symmetric, the identifiers are assumed to be ordered, i.e. there 
is only ddl223 but not dd3212. 
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